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Abstract: Starting from a commensurate triangular thin solid strip, confined within two hard struc- 
tureless walls, a stretch along its length introduces a rectangular distortion. Beyond a critical strain 
the solid fails through nucleation of "smectic"-like bands. We show using computer simulations 
and simple density functional based arguments, how a solid-smectic transition mediates the failure. 
Further, we show that the critical strain introducing failure is inversely proportional to the channel 
width i.e. thinner strips are stronger ! 
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I. INTRODUCTION 

Studies of small systems comprising of a few thousand 
molecules have become increasingly important with the 
advent of nano-technologyQ. In these mesoscopic length 
scales, a priori, there is no reason that mechanical behav- 
ior will be governed by continuum elasticity theory 0. In 
many situations one needs to study the effect of confine- 
ment on the structural and other properties of materials. 
Indeed, hard confinement, in one or more directions of- 
ten induce new interesting properties; for example, the 
rich phase behaviour of quasi two-dimensional colloidal 
solids 0, B confined between two glass plates showing 
square, triangular and "buckled" crystalline phases and 
re-entrant surface melting transitionlsl of colloidal hard 
spheres quite different from the bulk[^0,l3 behaviour. 

In an earlier paper |^ we have shown that a small 
confined solid in quasi one-dimension fails through nu- 
cleation of "smectic" like bands. This is very unlike 
a bulk solid, strained beyond it's critical limit, fail- 
ing through the nucleation and growth of cracks [ij, 

The interaction of dislocations or zones of plastic 
deformation Q] with the growing crack tip deter- 
mines the failure mechanism. Bulk solid show brittle 
or ductile|l|i failure depending on these interactions. On 
the other hand the two-dimensional confined solid shows 
ductile failure along with reversible plasticity — the plas- 
tic stress vanishes when the strain is removed and the 
fractured parts join up without any discontinuity. This 
is because of the high amount of orientational order im- 
posed by the confining walls which ensure that each smec- 
tic band is confined within a dislocation - antidislocation 
pair. On the removal of the strain these defects annihi- 
late in pairs leaving a perfect triangular lattice identical 
to the initial solid. 

In sectioninwe describe the phase structure of the two 
dimensional confined strip and possible phase transitions. 
We illustrate our discussion with previously published 9] 




FIG. 1: The confined solid is shown along with the centered 
rectangular (CR) unit cell. For an unstrained triangular lat- 
tice ax = ao and ay = \/3ao- Gi, G2 and G3 denote the 
directions of the three reciprocal lattice vectors (RLV). The 
third reciprocal lattice direction G3 is equivalent to the di- 
rection G2, even in presence of the walls. 

computer simulation data. In section ITTTl we present sim- 
ple density functional|l7l arguments showing the dif- 
ference between a bulk solid, confined-solid and a smectic 
phase as well as a Lindeman-like criterion for nucle- 
ation of the smectic. In section IIVI we give the results 
and discuss their implications. We end by offering some 
concluding remarks and discussions on future directions 
of research in section Ivl 



II. THE SYSTEM 

A two dimensional bulk solid made of molecules in- 
teracting via spherically symmetric potentials always 
stabilizes in the triangular lattice configuration. 
For specificity and simplicity we consider hard -disk 
molecules, which are restricted from overlapping with 
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each other due to an infinitely large energy cost and re- 
mains non-interacting when they do not overlap. Apart 
from being easily accessible to theoretical treatment j2r|. 
experiments with nearly "hard" (macro-) molecules viz. 
sterically stabilized colloids(21| are possible. The hard - 
disk free energy is entirely entropic in origin and the only 
thermodynamically relevant variable is the number den- 
sity p = N/V or the packing fraction 77 = (7r/4)pd^. The 
energy scale of the system is entirely set by fc^T where ks 
is the Boltzmann constant and T the (kinetic) tempera- 
ture. Accurate computer simulation 0, vindicated the 
defect unbinding theory 10] of melting of the hard -disk 
solid below -qm = .706. 

Imagine a narrow channel in two dimensions of width 
Ly defined by hard walls at y = and Ly {Vviai\{y) — 
for d/2 < y < Ly — d/2 and = 00 otherwise) and length 
Lx with Lx ^ Ly. Therefore, ni layers of an undistorted 
triangular solid of lattice parameter oq and diameter d 
may be accommodated (Fig^ if Ly is commensurate |^ 
with the solid packing fraction, i.e. 



— {ni - l)aQ + d. 



(1) 



For a system with a constant number of particles N 
and Ly, oq is decided by the packing fraction 77 alone. 
Note that L^ — rixtto — Nao/ni, and oo is given by 
p = N/LxLy. Defining x{v,Ly) = 1 + 2{Ly - d)/V3ao, 
the above condition reads x = integer = n/. Violation of 
Eqn. induces an internal rectangular strain measured 
from a reference triangular lattice of ni layers; 



The external strain, on the other hand, imposed by 
rescaling L^ keeping Ly fixed, e — (r]Q — ri)/r] is measured 
using the initial triangular lattice as the reference state. 
As the solid is free to move locally, it may readjust it- 
self to increase (or decrease) the number of layers ni in 
response to the external strain ensuring that the inter- 
nal deviatoric strain remains small (minimization of free 
energy). It can be shown that 9], 
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ni 
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(3) 



where the number of layers n/ is the nearest integer to 
X so that Ed has a discontinuity at half-integral values of 
X. This internal strain Ed is related non-linearly to e and 
may remain small even if e is large (Fig|21). 



(a) 



(b) 




FIG. 3: Structure factors of different phases, (a), (b) and (c) 
are the structure factors corresponding to solid, smectic and 
modulated liquid phases respectively. 
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The lattice parameters of the strained, centered - 
rectangular (CR) unit cell a.^ and ay are shown in Fig^ 
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FIG. 2: A plot of internal strain Sd as a function of external 
strain e. The jumps in Sd corresponds to half-integral values 
of X- 



Walls induce density modulations, with peaks and 
troughs running parallel to the walls. The ensemble av- 
erage of the local density p(r) =< n ~ ^'-^ > 
shows asymmetric (ellipsoidal) density profiles at the 
lattice points, with the semi-major axes lying along 
the walls. The structure factor is defined as pc = 



7 Yfj,k=i exp(-«G.rjfe)^ . For G = ±01(77), the re- 
ciprocal lattice vector (RLV) correspond to the set of 
close-packed lattice planes of the CR lattice perpendic- 
ular to the wall, and for ±62(7?) and ±63(77) the four 
equivalent RLVs for close-packed planes at an angle ( 
= 7r/3 and 27r/3 in the triangular lattice) to the wall 
(see Fig. It is useful to define the Lindemann pa- 



rameter I —< {u^ 



where the angular brackets denote averages over config- 
urations, i and j are nearest neighbors and u"i is the a- 
th component of the displacement of particle i from it's 
mean position. The parameter / diverges at the melting 
transition 1^. 

Apart from the solid phase (Fig^ and Fig|3Ia)), the 
externally imposed density modulation may give rise to 
a modulated liquid or a smectic phase. Structurally a 
smectic is made up of overlapping asymmetric local den- 
sity profiles along the walls generating continuous strips 
of density maxima running parallel to the walls. Peaks in 
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the structure factor corresponding to ±G2(??) vanish, al- 
though for G — ±61(77), one continues to obtain strong 
peaks Fig|2Ib). In this phase, however, inter-layer parti- 
cle exchanges are suppressed causing the Lindemann pa- 
rameter to remain small. In a modulated liquid, on the 
other hand, this interlayer exchange is large. The struc- 
ture factor of a modulated liquid (FiglSJc)) features a 
ring-like maximum characteristic of liquid, in addition to 
somewhat strong peaks corresponding to Gi. Transitions 
among these phases are observed as the external strain 
is imposed strating from a perfect triangular solid. The 
sequence of phase changes is shown in Fig01 the result 
of extensive Monte Carlo simulations of this system|9(|. 
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FIG. 4: Results of Monte Carlo (MC) simulations of = 650 
hard -disks confined between two hard walls separated by a 
distance Ly = 9 d. For each 77, the system was equilibrated 
over 10® MC steps (MCS) and data averaged over a further 
10*^ MCS. At 77 = 0.85 we have a strain free triangular lattice. 
Plots show the structure factors pG^^i ~ 2(o) for RLVs 

Gi(77), averaged over symmetry related directions, as a func- 
tion of 77 and the Lindemann parameter ?(□). The lines in the 
figure are a guide to the eye. At 77 « .78 pa^ jumps to zero 
while pd shows a small dip. The region with .7 < rj < .78, 
therefore, contains a smectic phase which melts for 77 < .7, at 
which point the Lindemann parameter / diverges. 



about by this overlap in the local density. Since me- 
chanical failure in our system is a consequence of a phase 
transition, it is reversible — as the strain is reduced back 
to zero, the stress also vanishes and the perfect triangular 
lattice is recovered 1^. 

Within density functional theory the excess grand 
potential of a non-uniform liquid containing a density 
modulation p{r) over the uniform liquid of density pi is 
given by. 



An 

1^ 



= J dr[p{r)log{p{r)/pi)-Sp{r)] 



I dv'C{\v~v'\)5p{v)5p{v'). (4) 



Here 5p(r) = p(r) — pi and C(r) is the direct corre- 
lation function of the uniform liauid[20j. A functional 
minimization of the free energy yeilds the following self- 
consistency equation for the density: 



P(£) 
Pi 



exp 



dr'C(|r-r'|)(5p(r')] 



(5) 



In principle one should solve the above equation within 
the constraints imposed by the walls and obtain the equi- 
librium p{y). Substitution of this p(r) into Eqn^l gives 
the equilibrium free energy and phase transitions. While 
we intend to carry out this procedure in the future, we 
must point out that for the present problem, this is com- 
plicated by surface terms and anisotropic, external, fields 
which are difficult to incorporate. In this paper we shall 
take a much simpler route in exploring the various con- 
ditions for the solid -smectic transition given the nature 
of the (the order parameters) obtained from our sim- 
ulations. 

One may expand, therefore, the logarithm of the local 
density profile logp(r) in a Fourier series around a 
lattice point at the origin, to get. 



III. THEORY 



p(r) = A/- exp [ 2Co ^ pg, cos(Gi.r) ) (6) 



As shown earlierj^, the failure of commensurate solid 
under tensile strain imposed in the manner discussed in 
the previous section, comes about through the nucleation 
of smectic bands within the solid. Monte-Carlo simula- 
tions show, at half-integral x where the local internal 
strain Sd becomes discontinuous, p(r) at nearest neigh- 
bour sites overlap along the x-direction, parallel to the 
walls, generating smectic bands. The stress associated 
with Ed vanishes at these points and the solid fails under 
tension. In this section we shall show, using simple den- 
sity functional^, arguments, that the phase transi- 
tion and the consequent tensile failure (a smectic cannot 
support stress parallel to the smectic layers) is brought 



where Cq is a constant, of order unity, denoting the 
Fourier transform of the direct correlation function cal- 
culated at a q-vector corresponding to the smallest RLV 
set of the solid. We have kept contributions only from 
this set. 

For a perfect triangular lattice, the RLVs are Gi = 



y'¥- sm( 



f) and G3 



— ao. Using these 



yf, = cos(f) , 

*fcos(f)-y|^sin(f), where - ^ 
relations and the fact that in the presence of confining 
walls the Fourier amplitudes denoting solid order are vir- 
tually constant upto the transition and = PG3 7^ 
PGi, EqnEl gives, 
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p(r) =AAexp{Co(2pGi +4pGj}cxp (-^Cq ('^^ {(2pGi + Pcjy' + SpG^a;^} 



(7) 



Clearly the density profile is Gaussian, of the form, 
p(r) ^ exp{—y'^ /2ay — x^/2al). Therefore, the spreads 
of density profile in x and y-directions are given by ax 
and ay respectively, with 
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(8) 
(9) 



In the absence of walls, pgi = PGq. niaking a^ = o'y, 
i.e. the density profile comes out to be symmetric in both 
directions, as expected for the bulk triangular solid. The 
presence of walls make ay < a^ making the density profile 
elliptical with larger spread in x-direction, the direction 
parallel to the walls. Two neighbouring density profiles 
will overlap to form a smectic if ax > cix- This leads us 
to the definition of a measure of overlap Oi = (ax/ax). 
The condition O; > 1 is then the Lindcmann criterion 
for nucleation of the smectic phase. Remembering Ux — 
ao{ni - - 1), we get. 



O, 



1 



1 
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^yCoPG^ m - 1 



(10) 



Whenever, pc,^ i.e. with the loss of solid order 
Oi diverges although ay remains finite, since pc^ ^ 
in presence of the walls. This indicates a solid-smectic 
transition. However, even before pc^ the quantity 
A — and therefore Oi shows large jumps at those 

internal strain (e^) values where x becomes half-integer. 
It is interesting to note that, at these points Sd has dis- 
continuities and the system failsQ- This shows that the 
mode of failure predicted by our theory is through a solid- 
smectic transition. The fact that pGi remains non-zero 
even at very small densities, due to the confinement from 
the walls, gives rise to the density modulations in the 
confined liquid. 

We have shown therefore that the overlap in the den- 
sity profiles may be used as an "order parameter" for the 
solid to smectic transition. We show below that jumps in 
this order parameter tantamounts to mechanical failure 
of the solid. 



IV. RESULTS AND DISCUSSION 

We begin with studying the overlap A as a function of 
external strain e. For specificity, we start from a triangu- 
lar solid of packing fraction rj = .85 with Ly commensu- 
rate with a n/ = 10 layered solid. With increasing strain 




FIG. 5: For a 10-layered solid with Ly commensurate with 
the initial strainless triangular structure at 77 = .85 overlap 
term A is plotted as a function of external strain e. Density 
profile overlap shows a jump increase at strains e > .1, the 
failure strain value0. 



initially the overlap A reduces due to increased separa- 
tion (ttx) between neighbouring lattice points. But above 
a strain (e) of about 10%, x reaches the half-integral mark 
and A shows a discontinuous increase, indicating large 
overlap between neighbouring density profiles along the 
wall ; indicating a solid to smectic transition (FiglSJ at 
the failure strain e*. With further increase in strain the 
overlap reduces, again due to increased separation be- 
tween neighbouring lattice points. At higher strains the 
smectic melts into a modulated liquid due to increased 
fluctuations connected with the reduced density 0. 

We have performed this calculation for various Ly 
values commensurate with starting triangular solids of 
ni = 2, 3... 20 layers at packing fraction 77 = .85. We 
found out the failure strains e* at each Ly and plot- 
ted them in FigEl as a function of Ly. This clearly 
shows that the failure strain reduces with increase in Ly. 
This demonstrates the fact, derived earlier from Monte- 
Carlo simulations Jij, that thinner (smaller Ly) strips are 
stronger! 

In Fig|3we have plotted the overlap term A with in- 
creasing interwall separation Ly at rj = .85. The jumps, 
as usual, indicate failure strains corresponding to discon- 
tinuities in the internal strain Ed at half-integral values 
of X- The plot shows that the amount of overlap at the 
failure strains e* reduces with increasing Ly indicating 
that at large interwall separations the system starts to 
behave as a bulk solid and more conventional modes of 
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FIG. 6: Failure strains e* for various interwall separations Ly 
confining n; = 2 — > 20 layered triangular strips at = .85 is 
plotted. Failure strain decreases with increase in Ly. 




failure viz. through formation and interaction of cracks 
and twin boundaries starts becoming active. 



V. CONCLUSION 

In this paper we have described how mechanical and 
phase behaviour are intimately related to each other for 
a narrow solid strip confined between hard walls. Both 
may therefore be described using simple arguments based 
on density functional theory. The particles of the solid 
are assumed to interact via hard -disk interactions Apart 
from constrained hard -sphere colloids plj where our re- 
sults are directly testable, a similar fracture mechanism 
may be observable in experiments on the deformation of 
mono-layer nano beams or strips of real materials pro- 
vided the confining channel is made of a material which 
is harder and has a much smaller atomic size than that 
of the strip m). 

We believe that phase transition in confined systems, 
as described in this paper, may have future applications 
in nano-technology. The feature of reversible plasticity^ 
may be used to produce nano scale mechanical switches. 
For example, thermal conductivity and electrical conduc- 
tivity of such narrow strips after nucleation of smectic 
bands are expected to reduce drastically. These proper- 
ties need to be investigated in detail. Due to reversible 
plasticity a small change in strain in opposite directions 
can vary these properties of confined nano-strips giving 
rise to a possible use of it as mechanical switch in nano- 
devices. 
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